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ABSTRACT 

We forecast constraints on primordial non-Gaussianity achievable from forthcoming surveys by 
exploiting the scale-dependent halo bias introduced on large scales by non-Gaussian initial conditions. 
We explore the performance of exploiting both the shape of the galaxy power-spectrum on large scales 
and the cross-correlation of galaxies with Cosmic Microwave Background maps through the Integrated 
Sachs Wolfe effect. We find that future surveys can detect primordial non-Gaussianity of the local 
form with a non-Gaussianity parameter |/nl| of order unity. This is particularly exciting because, 
while the simplest single-field slow- roll models of inflation predict a primordial |/nl| *C 1, this signal 
sources extra contributions to the effective /nl of large-scale structures that are expected to be above 
our predicted detection threshold. 

Subject headings: cosmology: theory, large-scale structure of universe - galaxies: clusters: general - 
galaxies: halos 
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1. INTRODUCTION 

Tests of deviation from Gaussian initial conditions of- 
fer an important window into the very early universe and 
a powerful test for the mechanism which generated pri- 
mordial perturbations. 

While standard single-field models of slow-roll inflation 
lead to small departures from Gaussianity, non-standard 
scen arios allow for a larg er level of non-Gaussianity (see, 
e.g. jBartolo et all (|2004D and references therein). 

The standard observables to constrain non-Gaussianity 
are the Cosmic Microwave Background (CMB) and the 
Large-Scale Structure (LSS) of the universe. Tradition- 
ally, the most popular method to detect primordial non- 
Gaussianity (NG) has been to measure the bispectrum 
or the three-point function of the CMB dVerde et all 
[20001: IKomatsu et all [20051 : lYadav fc Wandem 12008ft . as 
the LSS bispe ctrum is sensitive to primordial NG only at 
high redshift (IVerde et a l. 2000; Scoccim arro et al.ll2004t 
Sefusatti fc Komatsull2007l : ICoorayl 12006: Pillcpic h et al.l 
2007). 

A powerful technique is based on t he abundance 
(iMatarrese et alJl2000t IVerde et a,l.ll2001t lLoVerde et all 
120071: IRobinson fc Baker! l2000t IRobinson et all I2000D 
and clustering dGrinstein fc Wisdll986t IMatarrese et al.l 
1986| iLucchin et aIT !l988) of rare events such as dark 
matter density peaks, as they trace the tail of the 
underlying distribution. These theoretical predictions 
have been tested against numerical N-body simulation s 
dKang et al.ll2007t iGrossi et al J 120071: IDalal et alJ l2007D. 

iDalal et all (|2007|) : IMatarrese fc Verdel ||2008T) showed 
that primordial NG affects the clustering of dark mat- 
ter halos inducing a scale-dependent large-scale bias. 
Here we argue that this effect could be used to con- 
strain NG through the Integrated Sachs Wolfe (ISW) 
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effect ([Sachs fc Wolfd [1967) and through the shape of 
the galaxy power-spectrum on large scales. We fore- 
cast how future galaxy surveys could constrain Gaus- 
sianity via this "halo-bias" effect. We find that con- 
straints from surveys which provide a large sample of 
galaxies or galaxy clusters over a volume comparable to 
the horizon size (e.g., DES, PanSTARRS, LSST, EU- 
CLID, ADEPT) are competitive with CMB bispectrum 
constraints achievable with an id eal CMB experiment. In 
particular, iBartolo et al.l (|2005l ) showed that even when 
the primordial /nl is tiny, the evolution of perturbations 
on super-Hubble scales, yields extra contributions to the 
effective /nl relevant for the LSS, which are configura- 
tion and redshift dependent. These contributions are of 
amplitude comparable to the forecasted errors and can 
therefore no longer be neglected. 

Alon g the way we offer physic al insights in the find- 
ings of IMatarrese fc Verdel (|2008fl (hereafter MV08) and 
explai n the connections to the approach of IDalal et al.l 
(|200l : we show that only the formalism of MV08 can 
be used to correctly handle the non-linear contributions 
to the primordial NG. 

2. NON-GAUSSIAN HALO BIAS IN CONTEXT 



IMatarrese fc Verdd (|2008l ) generalized to NG ini- 
tial con ditions results and techniques developed in 
the 80s dGrinstein fc Wise! Il986t IMatarrese et all 119861 : 
ILucchin et al. 1988) to relate the clustering properties of 
the collapsed structures (halos) to those of the under- 
lying dark matter distribution for Gaussian initial con- 
ditions. This approach, which we briefly review here, 
yields an analytic expression for the bias of dark mat- 
ter halos for non-local and scale-dependent NG models 
in which the primordial bispectrum of the potential is 
the dominant higher-order correlation and has a gen- 
eral form. The starting point of MV08 is to consider 
the expression for the correlation function of regions 
above a high threshold in the general NG case, which 
has the form dGrinstein fc Wisd 119861 : IMatarrese et al.1 
fl9MILTicchin et al.lll988h~ ~~ 



ZhMr) = -l + exp[X(r)] 



(1) 
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where X is a complicated expression that depends on all 
the n-point correlations of the underlying density field 
filtered on the mass-scale M, as well as on the threshold 
height (see MV08 for details). Then, one interprets the 
region above high thresholds as halos, and, for large sepa- 
rations r (small values of X), expands the exponential to 
first order and considers only terms up to the three-point 
correlation function. Finally, one Fourier transforms 
to obtain an approximated expression for the power- 
spectrum of halos in NG models. For NG of the type 
(iSalopek fc Bondfl990t iGangui et al.l ll99l iVerdeetaTl 
I200nt iKomatsu fc Sperge]ll2001h 



$ = + /nl(</> 2 -(^ 2 }) , 



(2) 



where $ denotes Bardeen's gauge-invariant potential 5 
and <fi denotes a Gaussian random field, the halo power- 
spectrum has the form 



P h (k,z) 



S 2 (z)P ss (k lZ ) 



[1 + 4/ NL <5 c (z)a(fe)] . (3) 



Here a(k) is the quantity in Fig. 3 of MV08, which 
is cx 1/k 2 at large scales; cr^ is the mass variance lin- 
early extrapolated to z — 0; S c (z) — A c (z)/D(z), with 
A c (z) the linear overdensity for spherical collapse (which 
weakly depends on z in non-Einstein de Sitter cosmolo- 
gies), and D(z) linear growth- factor of density fluctua- 
tions normalized to unity at z = 0. 

The Lagrangian bias &z, of the halos (all correlations 
and peaks considered here are those of the initial den- 
sity field, linearly extrapolated till the present time) is 
then b Lth (z,M) ~ bf h (z,M)[l + 2f Nh 5 c (z)a(k)\, with 
h (z, M) ~ 5 c {z)/{a 2 M D{z)). Making the standard as- 
sumption that halos move coherently with the underly- 
ing dark matter, one obtains the final Eulerian bias as 
bs = 1 + 6l, so that 



= 1 



' D{ Zo ) 



[l + 2f Nh S c (z f )a(k)] , (4) 



M 



where, following the approach of I Cat elan et al.l {1998), 
we have made explicit the dependence on both the halo 
formation redshift Zf and the observation redshift z Q . 
This expression for the NG bias of halos of mass M 
is scale-dependent and increases rapidly at large scales. 
The approach of identifying peaks with halos is valid 
for rare (massive) halos. Indeed, the above expression 
for the Gaussian Lagrangian h alo bias is approx i mate ; 
a more accurate expression is ([Efstathiou et al. 1 119881 : 



a more accura te expression is (Mstatn 
iCole fc Kaiserlll989HMo fc Whitelll996f) 



b? h (z ,M,z f ) 



D(z ) 



Sc(zf) 



M 



5 c (zf) 



(5) 



In addition, and as discussed below, this expression can 
be further refined with the help of N-body simulations. 
For objects that did not undergo recent mergers, Zf 3> 
z , the bias is well approximated by Eq. (j4|). Eq. ([4} 
however also applies to the case Zf k, z$ (rapid mergers) 

5 On sub-Hubble scales, Bardeen's gauge-invariant potential <t> 
reduces to the usual Newtonian peculiar gravitational potential, up 
to a minus sign. In the literature, there are two conventions for 
Eq. (0: the LSS and the CMB one. Following the LSS convention, 
here <E> is linearly extrapolated at z = 0. In the CMB convention 
<I> is instead primordial: thus /nl = (g{z = oo)/g(0))f^ IB . 



for S 2 3> <J 2 M , i.e. large masses and/or high formation 
redshifts. 

Note that Eq. g]) can be rewritten as &{ N £ = bf h + Ab, 
where 

Ab = 2f NL S c (z f ) (bf{z , M, z f ) - 1) a(k) . (6) 

One may note that for /nl large and negative, Eq. ([3]) 
would formally yield &/ l NL and Ph{k) negative on large 
enough scales. This is a manifestation of the break- 
down of the approximations made: a) all correlations 
of higher order than the bispectrum were neglected: for 
large NG this truncation may not hold; b) The expo- 
nential in Eq. ([T]) was expanded to linear order. This 
however could be easily corrected for, remembering that 
the P(k) obtained in Eq. ((3|) is in reality the Fourier 
transform of X, the argument of the exponential. One 
would then compute the halo correlation function using 
Eq. {T]) and Fourier transforming back to obtain the halo 
p ower-spectrum. 

iDalal et all (|2007f ) and ISlosar et all (|2008f) use the 
peak-background formalism to o btain an equation s imilar 
to Eq. (j4|). In particular, while IDalal et al.1 (|2007f ) relies 
on the spherical col lapse model (and thus o n the stan- 
dard Press-Schecte r (|Press fc Schechtei1ll974j ) approach), 
ISlosar et al.l (|2008f l extend and reformulate it so that it 
relies on the extended Press-Schecter approach and on 
the universality of the mass function, and can be ob- 
tained for any mass function, even one that is a fit to 
N-body sim ulations. Th e adva ntage of their formulation 
over that of IDalal et all (|2007n is that they can include 
the description of the effect of halo mergers. The effect 
of halo merger can be analytically described only for the 
standard Press-Shechter mass function, but they argue 
that the scaling of the correction for mergers could be 
calibrated from N-body simulations. 

We note here that their derivation can only be carried 
out for NG of the local type. On the other hand, the for- 
mulation of MV08 is more general: the extended Press- 
Shechter approach, halo mergers and mass functions that 
are better fit to N-body simulations than the standard 
Press-Shechter can be readily taken into account by sub- 
stituting b'f h in Eq. ([5]) by the peak-background split 

bias which is obtained as b ( ^ h = —n~ 1 dn/d5 c , where n 
denotes the halo mass function for Gaus sian initial condi- 
tions, which could be given e.g. by the lSheth fc Tormenl 
( 1999) formula or by a fit to simulations. 
In particular, Eq. ([5]) is replaced by 



b% h (Za,M, Zf) = 



D(z ) 



q5 c {zj 



M 



S c (zf) 



(7) 



+ 



2p 



5 c (z f )D(z ) 



1 



M 



The parameters q and p account for non-spherical col- 
lapse an d fit to numerical simulations yield q ~ 0.75, 
p = 0.3 (|Sheth fc Tormenl 1 19991: ISheth et al.ll2001f ). The 
correction for non-spherical collapse also applies to the 
NG correction to the halo bias: 



A6 = 2q'f NL 6 c {b% - l)a(k) 



(8) 



where q' can be calib rated to N-body si mulations and is 
found to be q' w 0.8 (|Grossi et al.ll2008D . We stress here 
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that it is particularly important to be able to account 
for general non-local and scale-dependent NG character- 
ized by a given bispectrum of the potential. Indeed, as 
shown bv lBartolo et alj (|2005f) . there are extra contribu- 
tions to the bispectrum that come in at the same level 
as the primordial signal: in other words the primordial 
contribution is enhanced for LSS, yielding configuration 
and redshift dependent contributions to the effective /nl 
which cannot be neglected. Note that the NG bias for- 
mula of MV08 is fully general and can easily account for 
non-constant /nl- 



3. METHOD 

Before we present the forecasts we need to make two 
considerations. 

Halos vs Galaxies: The theory developed in MV08 
and above describes the clustering properties of halos, 
but we observe galaxies. Different galaxy populations 
occupy dark matter halos following dif ferent prescrip- 
tions . If we think in the halo - model (jCoorav fc Shethl 
2002t iScoccimarro et all 120011: iPeacock fc Smith! 120001 : 
Smith et alj 120031 ) framework, at very large scales (as 
those relevant for this analysis) , only the "two- halo" con- 
tribution matters and the details of the halo occupation 
distribution of galaxies (the so-called "one-halo" term) 
is unimportant. In particular, the galaxy population 
known as "luminous red galaxies" (LRG), is known to 
be old and free from recent merger activity. For this 
population, the modeling of MV08 should offer a good 
description. Emission lines galaxies on the other hand 
may be affected by recent merger activity and so their 
6^ may need to be modified, as described in §2. As we 
discuss below, with this modification, the effect of uncer- 
tainties in the Gaussian bias enters in our estimates only 
through the shot-noise contribution to the signal-to-noise 
ratio. 

Detection vs Measurement: Before we present the fore- 
casts we need to make an important distinction between 
"detection" and "measurement" . To compute the statis- 
tical significance of a detection we need to compute the 
significance of deviations from the null hypothesis: in 
particular, the fiducial model used in the calculation has 
/nl = and the error bars are also computed assuming 
the null hypothesis /nl = 0. However to carry out a 
measurement of /nl ^ 0, the theoretical model and the 
error bars must be computed as functions of /nl- Here 
we will report forecasts for detection of /nl- 

We consider two probes: the large-scale power- 
spectrum of galaxies and the ISW effect. In both cases, 
since we will consider tracers of rare halos, we will as- 
sume Zf = z , yielding a possibly conservative estimate 
of the errors on /nl- 



3.1. Forecasts from the shape of the large-scale 
power- spectrum 

In the Fisher matrix approach to error forecasts we can 
write 



ln£ 



1 (AP(fc)) 5 



1 (P(fc)4<5 c (z)a(fc)/ NL ) 2 



(9) 



where we have assumed a Gaussian likelihood 6 . An 
estimate of the error in /nl, <?7 NL (fc), at a given k is 
given by 

d 2 \\nC\ 16P 2 a(k) 2 S 2 , . 

— ^ — - (10) 



1 



The relative error in P for a shell in /c-space of width Afc 
and for a survey with effective volume V e ff is 

2 2 

(11) 



\ P J ATrk 2 AkV eff /(2TT 3 ) 



where V e ff = V(l + l/(nP)), with V the survey volume, 
and fi the average density of galaxies. If shot-noise is 
sub-dominant (i.e. hP » 1) V e fj = V. 
The total error from a fc-range from k m i n to k max is: 

8 A c (z) 2 



1 



> r 



2vr 2 D(z) 2 



■V, 



eff 



l{k) 2 k 2 dk. (12) 



If we divide the survey in redshift slices centered around 
Zi then the errors obtained combining different redshift 
slices (if uncorrelated) is 

1^1 

(13) 



' F 
JN 



We set k max to be 0.03 h/Mpc and k min to be greater 
than the where V is the volume of the shell 

considered (Az = 0.1). Note that the effect of NG 
alters the broad-band behavior of the P(k) on very 
large scales, which is not affected by the precision with 
which the radial positions of the galaxies in measure. 
Thus, we can treat photometric and spectroscopic sur- 
veys on the same footing. The requirement of surveying 
a large volume of the universe and sampling highly bi- 
ased galaxies to beat shot-noise, which is a key point 
for BAO surveys is also a bonus for constraining primor- 
dial NG. In particular, for the k m i n we use we find that 
Pg(k,nin) — Pg(k = 0.2/i/Mpc), thus the shot-noise re- 
quirement for BAO surveys of nP(k = 0.2h/Mpc) > 1 
implies that for all scales of interest here fiP 3> 1. We 
have checked that our results do not change if we impose 
nP ~ 3. 

While for BAO surveys measuring redshift accuratel y 
is crucial (jSeo fc Eisensteinl 120031 : Blake fc Bridldl2005l) . 
for this application is not important. 

3.2. Forecasts from the Integrated Sachs Wolfe effect 

The ISW effect probes the largest cosmological scales. 
As the NG effect goes oc 1/k 2 on large scales, this is a 
promising probe of NG. Here we follow lAfshordl ([2004) 
to quantify the significance of a detection of /nl through 
the estimate of the cross-correlation between the ISW 
effect with LRG galaxy distribution. For a galaxy sur- 
vey with the average comoving density distribution n c (r) 
as a function of the comoving distance r, in the Lim- 
ber approximation, the expected cross-correlation in the 
spherical harmonic space, can be written as 



C gT (i) = 



2 T 



J dr r 2 n c (r) 



dr n c (r)P&, g (k), (14) 



6 While strictly speaking the distribution of P{k) is non- 
Gaussian, this is a standard assumption in Fisher based ap- 
proaches. 
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TABLE 1 
Galaxy Surveys considered 



survey z range sq deg mean galaxy density (h/Mpc) 3 A/nl/<?' LSS 



SDSS LRG's 


0.16 < 


z < 0.47 


7.6 x 10 3 


1.36 


X 


10" 


-4 


40 


BOSS 


< ; 


! < 0.7 


10 4 


2.66 


x 


io- 


-4 


18 


WFMOS low z 


0.5 < 


z < 1.3 


2 x 10 3 


4.88 


x 


io- 


-4 


15 


WFMOS high z 


2.3 < 


z < 3.3 


3 x 10 2 


4.55 


x 


10" 


-4 


17 


ADEPT 


1 < 


z < 2 


2.8 x 10 4 


9.37 


x 


io- 


-4 


1.5 


EUCLID 


< 


z < 2 


2 x 10 4 


1.56 


X 


io- 


-3 


1.7 


DES 


0.2 < 


2 < 1.3 


5 x 10 3 


1.85 


x 


io- 


-3 


8 


PanSTARRS 


< z < 1.2 


3 x 10 4 


1.72 


X 


io- 


-3 


3.5 


LSST 


0.3 < 


z < 3.6 


3 x 10 4 


2.77 


X 


io- 


-3 


0.7 



where k — (£ + 1/2) /r, <f>' is the derivative of the gravita- 
tional potential with respect to the conformal time, and 
6g I £ Jn and T( m are the projected survey galaxy overden- 
sity and the CMB temperature in the spherical harmonic 
space, respectively. 

The expected dispersion in the cross-correlation signal 
is AC 2 gT (£) ~ C gg (£)C TT (£)[f sky (2£ + where / sky 

is the fraction of sky covered in the survey, and we as- 
sumed a small cross-correlation signal, i.e. C 2 T (£) <C 
C gg (£)C TT (£). 

For a galaxy distribution biased according to Eq. (|4|), 
dividing the survey in rcdshift shells, and following the 
same procedure of Section 3, the error in each shell at 
redshift z for a given I is 

_ 2 = -f[H(z)D(z)£((l + z)D(z))P ss (k,0)Ab(k,z)] 2 r 2 5r 
CT/NL " (21 + lfG TT (£)[P G (k, z) + n c (r)-!] 

(15) 

where 7 = 8/ sky (3TH 2 n m0 /c 3 ) 2 , k = (I + l/2)/r, 5r = 
{c/H(z))Az, Ab is Eq. © in the limit / NL = 1, and P G 
denotes the galaxy power spectrum in the Gaussian case. 
We impose fc min to be greater than the largest mode that 
can be sampled in each survey shell and k max = 0.03 
/i/Mpc. The total error is obtained summing up Eq. (|15p 
on all the multipoles I < 200 and integrating over the 
minimum and maximum redshift of each survey. 

For future large-scale galaxy surveys, we obtain 
A/nl = 7.6, 12.5, 11.5 for LSST, EUCLID and ADEPT, 
respectively. 



clustering approach to primordial NGis in principle more 
promising than the ISW one: the ISW signal is weighted 
at low redhift (z < 1), when dark energy dominates, while 
the effect of NG grows with redshift. However, the two 
approaches are affected by different systematics and thus 
should be considered complementary. 

It is interesting to compare the constraints on 
primordial NG achievable from the large-scale halo 
clustering with those achievable with the small-scale 
galaxy bispectrum. For example, comparing with 
ISefusatti &: Komatsul (|2007h we deduce that the halo- 
clustering constraints are a factor of 3 stronger than 
the bispectrum ones. The bispectrum however, through 
its dependence on the fc-space configuration, can be 
used to discriminate among different forms of NG. The 
CMB bispectrum for an id eal experiment ca n yield 
constraints of A/ NL = few (jYadav et alJ l2007h . The 
results of Table [1] indicate that constraints on /nl of 
order unity are achievable with future surveys, making 
it a highly competitive technique. We conclude that it 
is particularly important to be able to take into account 
general non-local and scale-dependent NG features 
characterized by a giv en bispectrum of th e potential. 
In fact, as shown by iBartolo et alJ (|2005l ). there are 
contributions to the bispectrum, which have a specific 
shape and redshift dependence and which come into 
play at the level of /nl ~ few. This is well above 
the detection threshold for forthcoming and proposed 
surveys, thus opening up the possibility to measure 
these secondary contributions to /nl- 



4. RESULTS & DISCUSSION 

Here we present forecasts of /nl constraints for forth- 
coming and future surveys. The surveys we consider and 
their specifications are reported in Table Q] along with 
the 1 — er error on /nl from the shape of the galaxy 
power-spectrum. The reported errors on /nl have been 
normalized by the correction factor for non-spherical col- 
lapse q' « 0.8. Note that the number of galaxies and 
the Gaussian bias enter in this signal-to-noise calcula- 
tion only through the contribution to the error due to 
shot-noise. The reported numbers are not dominated by 
shot-noise. 

This signal-to-noise calculation indicates that the halo 



While this work was being completed we became 
aware of Afsh ordi k, Tolley (|arXiv:0806.1061|) and of 
McDonald (arXiv:0806.1046). Our results are in good 
agreement with theirs. 
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